Critical fluctuations of time-dependent magnetization in a random-field Ising model 
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Cooperative behaviors near the disorder-induced critical point in a random field Ising model are 
numerically investigated by analyzing time-dependent magnetization in ordering processes from a 
special initial condition. We find that the intensity of fluctuations of time-dependent magnetization, 
x{t), attains a maximum value at a time t = r in a normal phase and that x( r ) an d t exhibit 
divergences near the disorder-induced critical point. Furthermore, spin configurations around the 
time t are characterized by a length scale, which also exhibits a divergence near the critical point. 
We estimate the critical exponents that characterize these power-law divergences by using a finite- 
size scaling method. 
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I. INTRODUCTION 

It has been known that earthquakes [l[ , acoustic emis- 
sions in a deformed complex material 2JL and Barkhausen 
noise in random magnets 0, [H, @, 0> HiEl exhibit distinc- 
tive power-law behaviors. All these systems possess some 
disorder and also they are driven by a slowly varying ex- 
ternal field. As a simple model of such systems, a ran- 
dom field Ising model (RFIM) under a slowly varying 
magnetic field has been investigated in order to elucidate 
the essential mechanism of the power-law behaviors. By 
performing numerical experiments of the RFIM under a 
slowly increasing magnetic field, it was found that the 
size distribution of avalanches, each of which represents 
a spin flipping in a connected region, becomes a power- 
law function at the critical strength of the disorder [J]. 
Then, the critical magnetic field at which the avalanche 
size becomes a system size is called disorder-induced crit- 
ical point. Since this phenomenon occurs due to the ex- 
istence of disorder, such power-law behaviors are called 
disorder-induced critical phenomena. 

Despite the extensive studies for the critical phenom- 
ena, fluctuations of magnetization, which might be the 
most naive quantity characterizing the criticality, has 
never been investigated. Related to this issue, it has 
been known that the magnetization as a function of 
the magnetic field shows the almost discontinuous be- 
havior at the disorder-induced critical point, although 
the precise determination of the transition type is still 
a controversy EH El- From this observation and 
considering the knowledge of conventional critical phe- 
nomena, static fluctuations of magnetization hardly ex- 
hibit singular behaviors near the disorder-induced critical 
point. These raise a naive question whether the disorder- 
induced critical point can be characterized in terms of 



fluctuations of time-dependent magnetization. 

In this paper, we present a positive answer to this ques- 
tion. Our key idea is to notice the recent extensive studies 
for critical fluctuations near an ergodicity breaking tran- 
sition in glassy systems QJ, El El El El El El • Here, 
we review these studies briefly. As an example, let us 
consider cooperative behaviors near the transition point 
in colloidal suspensions. In this system, near the tran- 
sition point, there exist long-range spatial correlations 
among movable particles during a time interval r, which 
is chosen as a typical relaxation time. When we intro- 
duce an appropriate quantity Q(r,t) that indicates the 
occurrence of large particle displacement at the position 
r during the time interval i, a discontinuous jump occurs 
in Q(r, t — > oo) accompanied with critical fluctuations of 
Q(r,r) at the ergodicity breaking transition. It should 
be noted that Q(r,t) characterizes the dynamical event 
intrinsic to glassy systems. 

These results motivate us to study the cooperative be- 
haviors near the disorder-induced critical point of the 
RFIM from the viewpoint of fluctuations of some dy- 
namical events. In particular, we consider ordering pro- 
cesses of the magnetization from a special initial condi- 
tion. We have numerically found that fluctuations of the 
time-dependent magnetization exhibit a critically diver- 
gent behavior near the disorder-induced critical point. 

This paper is organized as follows. In Sec. [Til we 
introduce the random-field Ising model and demonstrate 
its basic behaviors numerically. In Sec. IIII1 we address 
the main result of our study. Concretely, we present new 
critical exponents that characterize the phenomena near 
the disorder-induced critical point. The final section is 
devoted to concluding remarks. 
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Let us consider a three-dimensional cubic lattice A = 
{i = [x, y,z)\l < x,y,z < L}. A spin variable Oi G 
{ — 1, 1} is defined at each site is A. We study the RFIM 



2 



described by the Hamiltonian 

H = -J^2 GiGj - ^2(hi + h)ui, (1) 

where represents a nearest-neighbor pair of sites, h 
is a constant external field and the random field hi obeys 
a Gaussian distribution 



P(hi) = 



1 



V2vri? 2 



exp(- 



R_ 

2R 2 



(2) 



The time evolution of the spin variables is described by 
the following rule. We first choose a site at random. If 
the spin flip on the site makes the system energy lower, 
the sign of the spin variable is changed; otherwise the spin 
flip is rejected. At the next step, we choose a site at ran- 
dom again, and repeat the above-mentioned procedure. 
Here, a unit time is given by L 3 steps, which is called 
a Monte Carlo step per site (MCS). This time evolution 
rule is similar to traditional avalanche dynamics such as 
the one used in Ref. [6(, but our rule may be closer to 
standard Glauber dynamics. Indeed, it corresponds to 
the Metropolis method with the zero temperature. We 
should be careful to choose a time evolution rule particu- 
larly when we discuss finite temperature cases, but we do 
not enter a difficult question which time evolution rule is 
more physical. 

Here, we review the disorder-induced critical phenom- 
ena in the RFIM. Let us consider the quasi static change 
of the external field from h = — oo to h = oo. In this 
operation with a given R, h c (R) is defined as the spe- 
cial value of h (if it exists) at which a spin flipping in 
a region over the whole system first occurs. Then, R c 
is defined as the maximum value below which h c (R) ex- 
ists. Their actual values were numerically determined as 
(R c ,h c {R c )) ~ (2.16,1.44) @, which is called disorder- 
induced critical point. It was reported that power-law 
behaviors of the size distribution of avalanches were ob- 
served in the quasi static change of the magnetic field 
when R = R c . Such power-law behaviors are called the 
disorder-induced critical phenomena. 

In this paper, we study ordering processes of magneti- 
zation 



m(i) 



(3) 



from the initial condition m(0) = —1 in which all the 
spins are downward. We also focus on the case T — 
0, except for a brief discussion in the final section. As 
preliminary calculations, we measure 



m(i) = (rh(t)) 



(4) 



for several values of (R, h) , where (A) represents the aver- 
age of a physical quantity A with respect to the stochastic 
time evolution and quenched disorder. In the argument 
given below, we consider at least 20 samples of time evo- 
lution for each set of {hi} and at least 20 samples of {hi} 
for calculating average values. 
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FIG. 1: Ordering process of magnetization, m(t). L — 40 and 
R = 2.16 




FIG. 2: m(too) as a function of h for systems with three 
different sizes. 



Typical samples of m(t) for the case R = 2.16 are 
displayed in Fig. [TJ It is observed that the magnetiza- 
tion m(t) quickly approaches the equilibrium value m eq , 
which is approximately 1, from -1 when h is sufficiently 
large and that the ordering process becomes slower for 
the system with smaller h. Note that the equilibrium 
value m oq is slightly less than 1 because hi on some sites, 
whose values are largely negative, prevent the magnetiza- 
tion from approaching 1. Then, by decreasing h further, 
we find that there is a value below which m(t) does not 
reach the state with the equilibrium value m cq . 

In order to quantify this transition, we consider the 
quantity m(too), where is chosen as 1000 MCSs in 
our numerical experiments. We confirmed that the re- 
sults reported below did not depend on the choice of too 
when too > 1000 MCSs. (The lowest value of too depends 
on the system size.) On the basis of this, we define the 
frozen phase as the state with m(too) < mcq- Since a dis- 
continuous transition from m eq to m(too) is expected to 
occur at a certain value of h in the large size limit, we can 
replace the exact definition of the frozen phase with an 
operational one expressed as m(t oo) < 1— e, where e = 0.1 
for numerical simplicity. Indeed, as shown in Fig. [21 the 
transition becomes sharper when L is larger. Therefore 
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FIG. 3: Phase diagram. The circles represent frozen states. 
L = 40. 



FIG. 4: x(t) for three values of ft. R = 2.16. 



we expect that our operational definition provides an ac- 
curate determination of the phases in the large size limit. 
Presently, we fix L = 40 and test whether the system ex- 
hibits a frozen phase for several values of (R, ft) . The re- 
sult is summarized in Fig. [31 which suggests the presence 
of a transition curve hb(R, L, e) between the frozen and 
unfrozen phases. Note that the disorder-induced critical 
point (i?, c , ft c ) reported in Ref. @ appears to be located 
on this curve. More precise correspondences with previ- 
ous phase diagrams reported in Refs. [!, 0j| H3] will be 
discussed elsewhere. 



III. RESULT 

Now, we focus on the behaviors near the disorder- 
induced critical point (R c ,h c ). That is, by fixing R as 
i? c , we investigate the system with several values of ft 
near ft c . Following our motivation, we are interested in 
studying fluctuations of m(t). It should be noted that our 
study is concerned with fluctuations of relaxation events. 
Although such types of fluctuations of dynamical events 
have been studied extensively in an ergodicity breaking 
transition in glassy systems, to our knowledge, there have 
been no such arguments on statistical properties near the 
disorder- induced critical point. 

Since the simplest quantity characterizing the fluctua- 
tions of m(t) is given by 



X (t) = N (m(t) 2 ) - (m{t)Y 



(5) 



we first demonstrate the graphs of x(t) for a few values 
of ft in Fig. 2J It is observed that X (t) has a peak at a 
time t and that both Xm — x( T ) an d T increase when 
ft approaches ft c (— 1.44). Based on this observation, we 
next attempt to extract divergent behaviors of Xm and t 
by using a finite-size scaling analysis. 

Thus for the systems with L = 10, 20, and 40, we 
measured r(ft, L) as a function of (ft — ft c )/ft c . As shown 
in the inset of Fig. [5J when we plot r(ft, L)L~ a as a 



function of L x l e {h — h c )/h c , these three graphs do not 
depend on L, where a and 9 are fitting parameters whose 
values (0 ~ 0.7 and a ~ 1.7) are determined in such a 
manner that the three graphs are collapsed into a single 
curve as exactly as possible. Based on this result, we 
conjecture a scaling form 



r(ft, L) = L a F T 



ft — ft r 



L" 



(6) 



by using the scaling function F T . Considering the asymp- 



totic law F T (z) 



(with Q ~ 1.3) in the regime z 1, 



we expect the following critical behavior in the large size 
limit: 



(ft-ftc)" C - 



(7) 



In a manner similar to that in the analysis of r(ft, L), 
we assume a form of the finite-size scaling as follows: 



Xm (h,L) = L b F x 



ft - ft r 



(8) 



Indeed, from Fig. [5l we determine b and the scaling func- 
3.0 and we find the asymptotic rela- 



tion F x , where b 



tion F x (z) ~ z 1 (with 7 ~ 2.1) in the regime z 3> 1. 
We thus obtain 



Xm ^ (ft - ftc)" 



(9) 



These divergent behaviors observed near the disorder- 
induced critical point lead us to expect that some spins 
flip cooperatively around the time r. In order to describe 
the nature of the cooperative phenomena, we attempt to 
define a spatial correlation length that characterizes it, as 
carried out in studies on traditional critical phenomena. 
Here, it should be noted that the magnetization grows 
around the time r. In this paper, for numerical simplic- 
ity, we focus on the spin configurations at the time to 
such that m(to) = for each sample because the mag- 
netization is expected to grow at this time. We then 
consider the spatial pattern indicating whether the spin 
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FIG. 5: Finite-size scaling for Xm(h,L) and r(h,L) (inset). 
The three graphs for systems with different sizes are collapsed 
into a single curve. The statistical error bars are within the 
size of symbols. h c = 1.44 and R = 2.16. 
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FIG. 6: S(k)/S(0) as a function of k^{h,L) Inset: Finite- 
size scaling for £(h, L). Filled symbols represent those of the 
system for L=64 and the others for L=32. R — 2.16. 



on each site has already flipped by the time to. Thus, we 
measure 



p(j)=6(a j (t ),l), 



(10) 



where 6(m,n) represents Kronecker's delta. Calculating 
its Fourier transform 



K k ) = ^2 P(j) ■ j), 

jeA 



we define the structure function 



S(k)^-(\p(k)\ 2 ) 



(11) 



(12) 



where we set k = (fc,0, 0) in the argument given below. 
For S(k) obtained for several values of (h,L), we find 
that the fitting 




FIG. 7: Schematic phase diagram in the (R,H,T) space. 



is obtained well with n — 3/2, as shown in Fig. [6] This 
is called the Ornstein-Zernike form if n = 2 [12|]. Then, £ 
represents the correlation length characterizing the spa- 
tial pattern p(i). 

Now, in the same manner as those for r and Xm: w e 
perform a finite-size scaling analysis assuming the form 



H{h,L) = LF^ 



h — h c 



(14) 



The inset of Fig. [S] illustrates that this assumption is 
reasonable and that Fc(z) obeys the asymptotic relation 
F^(z) ~ z~ v for large z (y ~ 0.7). On the basis of this, we 
obtain the critical behavior of the dynamical correlation 
length: 



(15) 



Note that the obtained value v ~ 0.7 is close to the 
value 9 ~ 0.7, where is the exponent that character- 
izes the length scale appearing in the finite-size scaling 
method. This coincidence implies that the manner of the 
divergence for length scales is characterized by a single 
exponent. It should be noted that the exponents char- 
acterizing the length scale appearing in finite-size scaling 
analysis for statistical quantities related to avalanches 
were obtained in Refs. [13, HH , where the values of the 
exponent are close to that of 9 in our study. The relation 
among these results will be studied in the future. 

One may be afraid that another time evolution rule 
provides a different result, in particular with regard to 
the exponent £. Until now, we do not understand its de- 
pendency, but we conjecture that the value of v is not so 
influenced by the choice of the time evolution rule, be- 
cause of the consistency with the previous studies. These 
will be also studied in the future. 



IV. CONCLUDING REMARKS 



(13) 



We have presented the power-law divergences (JTJ) , ([H]) , 
and (|15p near the disorder-induced critical point in the 
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RFIM. We have determined the exponents by using the 
finite-size scaling method as 7 ~ 2.1, £ ~ 1.3 and v ~ 0.7. 
These values are regarded as preliminary values and more 
precise values will be determined by using systems with 
considerably larger sizes. At present, the accuracy of 
values is not of primary interest, but the existence of 
divergent fluctuations of time-dependent magnetization 
is rather important. Indeed, based on our results, we 
conclude that the phenomena near the disorder-induced 
critical point can be captured from the viewpoint of fluc- 
tuations of dynamical events. In this sense, the phenom- 
ena under consideration have common features with co- 
operative behavior in glassy and jamming systems. It is 
an important future subject whether the values of criti- 
cal exponents observed near the disorder-induced critical 
point are related to those in glassy and jamming systems. 

Last, we introduce a few examples of studies motivated 
by this conclusion. The first example is the phase dia- 
gram in the (R, H, T) space. Since the precise definition 
of the frozen phase appears to be complicated in the finite 
temperature case, we just plotted the region m(t oc ) < 
as a tentative frozen phase, using a Metropolis method. 
Figure [7] shows a schematic diagram of this phase. From 
this figure, one may recall the phase diagram of the jam- 
ming transition in granular systems (see Fig. 1 in Ref. 
[2^ | or Fig. 4 in Ref [Hj]). We also expect that a similar 



type of phase diagram can be obtained when we employ 
other time evolution rules such as a heat bath method. 
Thus this resemblance motivates us to study the com- 
mon aspects between granular systems and the present 
system. 

The second example is related to a theoretical frame- 
work. In addition to extensive analysis on the power-law 
distribution of avalanches 0,0, 9 , it was conjectured that 
the critical behaviors of avalanches in some s pin models 
with disorder are related to metastable states |24l. [25l [26| . 
It is interesting to investigate x(t) m such systems. Fur- 
thermore, we are interested in conducting a theoretical 
analysis of our numerical results. Since x(t) proposed in 
this paper has never been studied in the RFIM, such a 
theoretical analysis would shed light on a new aspect of 
the RFIM. 
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